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Abstract 

o 

<N We study the effect of resonances near the threshold of low energy (e) reactive scattering processes, and find an anomalous behavior 
of the i-wave cross sections. For reaction and inelastic processes, the cross section exhibits the energy dependence cr ~ e -3 ' 2 instead 
c3 of the standard Wigner’s law threshold behavior cr ~ e l,/2 . Wigner’s law is still valid as s —> 0, but in a narrow range of energies. 
yT We illustrate these effects with two reactive systems, a low-reactive system (FF + Cl) and a more reactive one (FF + F). We 
provide analytical expressions, and explain this anomalous behavior using the properties of the Jost functions. We also discuss the 
OI implication of the reaction rate coefficients behaving as K ~ l/T at low temperatures, instead of the expected constant rate of the 
Wigner regime in ultracold physics and chemistry. 
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G 1. Introduction 

<L> 

Ultracold gases allow a high level of control over the inter- 
^ action by using Feshbach resonances £lj] or by orienting ultra- 
^ cold molecules 01 . In addition to the study of phenomena in 
. degenerate quantum gases ranging from the BEC-BCS cross- 
^ over regime to solitons and multi-component condensates [4, 
i*|5t], such control also permits the investigation of exotic three- 
jf |. In i 


£^body Efimov states 
l——1 ing cold molecules 


arallel, the rapid advances in form- 
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„ „ nil have made possible studies 

of cold chemical reactions HlOl and their control fH . with 
application to a growing range of fields, such as quantum in¬ 
formation processing [ 1S lf5l llal ■ A key ingredient for these 
investigations are resonances near the scattering threshold. Al¬ 
though such resonances have been theoretically d and exper¬ 
imentally 118. 00 explored in the chemistry of low temper¬ 
ature systems, their effect on reaction rates is however not fully 
taken into account. 

In this article, we investigate how reaction and inelastic 
processes are affected by near threshold resonances (NTR) in 
the entrance channel of a reactive scattering system. We focus 
our attention on two benchmark systems containing hydrogen, 
namely FF+C1 and FF+F. We note that alkali hydrides 12111 . be¬ 
cause of their small mass, would also be interesting systems to 
investigate. By varying the mass of H, explore the modifica¬ 
tion of the behavior of the reaction cross section cr react from the 
Wigner’s threshold law cF 630 ' oc g -1 / 2 
ergy e. We generalize our initial treatment 
resonances in low enere 



clear 


and atomic |2f 


at ultralow en- 
of the effect of 
' collisions originally analyzed in nu- 
I collisions, to include a multi-channel 


formalism using Jost functions. 
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In Section [2) we review the general theory of reactive scat¬ 
tering and discuss the two benchmark systems FF+C1 and FF+F, 
and present the results of our calculations in Section [3] l or both 
systems. In Section [4] we give an explanation of those results 
based on the S-matrix, followed by an alternative description 
based on Jost functions in Section 0 Finally, we describe the 
equivalence of both approaches in Section [6] and give a simple 
physical picture in Section|7J before concluding in Section[8] 


2. Reactive scattering at low energies 

In this section, we review briefly the theory of low energy 
reactive scattering, and describe the two benchmark systems 
that we use to study the effect of s-wave near threshold reso¬ 
nances. 


2.1. General theory 


The scattering cross section from an initial internal state i to 
a final state / is given by [27] 




k 2 . 

i j =o 




’//I 


(1) 


e c 


where £ = \J - j \,...,/ + j and l' = \J - /1,...,/ + /; 
J-j + t-j'+t' is the total angular momentum, with molecu¬ 
lar rotational momentum j and orbital angular t in the entrance 
channel;, and corresponding quantum numbers J, j, and £ (the 
primes indicate the exit channel /). Here, e,- = /r/c 2 /2/j is the 
kinetic energy with respect to the entrance channel threshold, k, 
the wave number, and /./ = + m 2 the reduced mass in the 

entrance arrangement (with A standing for Cl or F). We are fo¬ 
cusing on the effect of resonances at ultralow temperatures, and 
so consider only s-wave scattering with C — 0, which requires 
J - j and thus (2/ + l)/(2 j + 1) = 1. In addition, we limit 
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ourselves to molecules initially in their rotational ground state 
(j = 0), such that we only need consider J - 0; consequently, 
the sum YiC i n 0} reduces to a single term: C = f. Hence, 
equation ([TJ simplifies significantly, and it now reads: 

cr f ^s i ) = ^\d fi -S J f f(k i )\ 2 . (2) 

i 

It is well known l28f l that in the zero-energy limit the cross 
sections can be expressed in terms of the complex scattering 
length aj - at - ifij, 

crj eact (e) = ^ ct / m (£) -> 

f*i 

of ast (e) = -> 


An —, 
k 

4n(a}+/%). 


(3) 


We obtain energy dependent rate coefficients by multiplying the 
cross sections with the relative velocity v re i = hk/ja. 



which can be thermally averaged over a Maxwellian velocity 
distribution to yield the true rate constant K(T). Note that the 
sum 2in Eqs. (j3j and 0 includes all possible final states, 
i.e., both reaction channels to form the products and quenching 
(inelastic) channels for the reactant. Thus, cr react and A' reacl de¬ 
fined above include all non-elastic outcomes of the scattering 
process, i.e., reaction plus quenching. We shall simplify our 
notation by omitting the channel subscript i (for s,-, kj, a /?,, 
etc.) and we also omit the superscript J in 5f. throughout the 
remainder of this article. 
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Figure 1: Rovibrational energy levels for FF and HC1. For H 2 we show the 
rotational series for v = 0 (dashed blue) and v = 1 (solid black). Similarly 
for HC1, we have v = 0 (dashed green), v = 1 (solid red) and v = 2 (dot- 
dashed cyan). The inset shows a sketch of the diagonal matrix element of the 
potential energy surface for the (v = \,j = 0) ro-vibrational entrance channel; 
the dashed line represents a weakly quasi-bound level for the van der Waals 
complex Cl- • • H 2 in the entrance channel H 20 ,0). 


We remark that, although the limits in Eq. (0 remain valid 
even in the presence of NTRs, their applicability will be limited 
to a much narrower domain of energies, and a new behavior 
will emerge for the remainder of the low energy regime, as we 
will show in this article. 


2.2. Systems with reaction barrier 

For our study, we selected two barrier-dominated reactive 
systems, namely Hi + Cl —> H + HC1 and Ht + F —» H + HF. 
Although both have reaction barriers, they have also different 
level structure energetics. 

As shown in Fig. [j] low energy scattering is purely elas¬ 
tic for Cl + H 2 with H 2 in its rovibrational ground state. We 
therefore consider the next vibrational level as an initial state, 
namely (v = 1, j = 0), such that the reaction to form HC1 be¬ 
come possible even in the ultracold regime. Note that quench¬ 
ing of H 2 by Cl will also be a possible (inelastic) outcome. 
Fig. [j] also shows the shallow potential well of the Cl- • • H 2 
van der Waals complex which supports quasibound states lead¬ 
ing to resonances in the entrance arrangement I 29 I 1 . This system 
was recently investigated at ultralow temperatures by Balakr- 
ishnan M and by us ll24ll . In this work, we use the potential 
energy surface developed by Bian and Werner & 

In the case of H 2 + F, reaction channels are already open, 
and the formation of HF at ultralow energies is possible even if 
H 2 is in its rovibrational ground state (v = 0, j = 0). Note that 
quenching is no longer possible in this case. However, there is a 
wide range of rovibrational states available for the product HF, 
as illustrated in Fig. [2] Again, the shallow potential well in the 
entrance channel can support quasibound states of the F- • • H 2 
van der Waals complex, and if such a state is near the threshold, 
it will produce resonance effects for low energy scattering. This 
system has been studied in the ultracold regime by Dalgarno 
and co-workers lEH. l33ll. We use the potential energy surface 
developed by Stark and Werner Q- 

As we did in our previous study of H 2 +CI [24J], we explore 
the effect of near threshold resonances (NTR) on reactive scat¬ 
tering by varying the mass of H. This was also done in |33!]; the 
channel thresholds in both arrangements H 2 -CI and H-HC1 (or 
H 2 -F and H-HF) shift at different rate, an approach similar to 
modifying the potential surface itself 13511 . 
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Figure 2: Same as Fig. [7]for HF. Only v = 0, j = 0, 2, 4 are shown for H 2 . 
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3. Results and discussion 


The results we present here were obtained using the abc 
reactive scattering code of Manolopoulos and coworkers l ITh ll, 
which we have optimized for ultralow energies in previous stud¬ 
ies of H 2 +D and H 2 +C1 Q. 


For the sake of completeness, the corresponding rate con¬ 
stants are plotted in Figs. [6ja) and (b). As m nears the reso¬ 
nance, AT react is significantly enhanced, scaling as T _i until T is 
small enough that the Wigner regime is reached and AT eact be¬ 
comes constant. The corresponding scaling for K tlasl changes 
from T~ l/2 for NTR to T i/2 for Wigner’s regime. 



Mass of H (u) 


3.1. Results for H 2 + Cl —> HCl + H 

Figure[3]shows the real and imaginary parts of the scattering 
length for H 2 (v = 1 ,j = 0)+Cl as a function of the mass m of 

H, with the open circles indicating the true masses of H and D. 
We note that H is located on the wing of a resonance, visible as 
a sharp increase for fi. 

The top panel of Fig. [4] shows the reaction cross section 
as a function of m for a few energies (in Kelvin). The simple 
cr oc oc e 1/2 scaling in Eq. ([3} implies equidistant curves 
for the energies chosen in the logarithmic scale. We find this to 
be true except near the resonances. This is more clearly illus¬ 
trated in the lower panel of Fig. [4] in which the rate constant is 
plotted; from Eq.©. we expect curves to coincide for all those 
four low energies, which is the case away from the resonances. 
However, as we near a resonance, the difference between the 
curves become more pronounced. 

In order to analyze this behavior, we selected three masses 
near the resonance, namely m = 1.0078 u= mu (true mass), 

I. 038 u, and 1.042 u; they appear as dashed vertical lines in the 
inset of the top panel of Fig. [4] We focus our attention on these 
three masses, and we analyze in detail the energy dependence 
of the cross sections. Fig. 0 a) shows the low energy behavior 
of cr ieact : as e —> 0, cr react reaches the Wigner regime, scaling 
as e 1 for all three masses, but for masses closer to the reso¬ 
nance, the scaling changes to e :,/2 . Fig. 0b) shows the elastic 
cross sections cr elast for the same masses; the Wigner regime’s 
constant cross section as s —> 0 changes to the expected e 1 
scaling for m near a resonance. 


Figure 4: Top panel: Reaction cross section for H 2 +CI as a function of mass m 
of H for different scattering energies in units of Kelvin. The inset shows cr react in 
the vicinity of the resonance with the vertical dashed lines corresponding to m = 
1.0078 u = mu (true mass, in black), 1.038 u (red), and 1.042 u (blue), from left 
to right respectively. Bottom panel: Rate constant for the same energies. 


3.2. Results for H 2 +F^>HF + H 

The results for the F + H 2 (0,0) reaction contain resonant 
features that are similar to those described above for the Cl + 
Hi (1,0) system. Figure [7] shows the reaction rate constant at 
T — 0, which is directly related to [3. see Eq. i|4ji. Although 
Fig- Eland Fig. [3] are similar, only two pronounced resonances 
are found for Hi + F, even though the mass of H varies over a 
wider range. Note that the real mass of H is again on the wing 
of a resonance. Dalgarno and co-workers [33] have obtained 
similar results, but their work has explored a limited mass range 
(only up to m = 1.5 u, with one additional data point for m = 
mo ~ 2 u) and they only found one resonance (near m = 1.1 u). 

The reaction cross section and rate constant follow the same 
behavior as the mass gets closer to the resonance. This is illus¬ 
trated in Fig. 0for cr react , with the NTR scaling (s~ 3/2 ) becom¬ 
ing more apparent as m gets closer to the resonance, while the 
Wigner regime (with e 1/2 scaling) is pushed to lower energies. 
The corresponding behavior for the rate constant is shown in 

Fig. m 

To explain the appearance of the NTR regime, we turn our 
attention to the analytical properties of the S -matrix. 


3 



























m ini 11 mini iiiiiml iiiiiml iiiiiml i.iimi i m i.ii iiiiim. iilmil iifliiu 

io " 9 io " 6 io " 3 10° 


Kinetic energy (K) 

Figure 5: Energy dependence of the reaction cross section (a) and elastic cross 
section (b). The results shown are for the entrance channel H 2 (v=l, y=0) + Cl 
—» H + HC1, and they correspond to different masses of H, as indicated for each 
curve. 


4. S-matrix approach for t = 0 

The anomalous behavior of the cross sections shown above— 
an abrupt increase followed by a gradual transition into the 
Wigner regime, as the energy decreases—is due to the presence 
of a resonance pole near the threshold of the entrance channel. 
To understand this, we pay attention to the S matrix. Recalling 
the unitarity condition, 2/ |5/,| 2 = |S,-,j 2 + Tif# 1 5 /,j 2 = 1, and 
using Eq. ©, we rewrite the cross sections in Eq. (0 explicitly 


in terms of the diagonal element S a(k): 


rr rcact (k) = ^(l _|.S’ 11 (k)| 2 ), 

(5) 

cr elast (k) = L\l - S tt (k )| 2 . 

(6) 


In this section we use the single channel case heuristically 
to guide our analysis of S „(&), while in the next section we will 
present a more rigorous approach based on Jost matrices for the 
general case of a many channel problem. 


4.1. Pole factorization; single channel case 

For a single channel problem, we express the S matrix for a 
given partial wave in terms of the Jost function Tr'- 


S { (k) = 


Tt(-k) 

m) ■ 


(7) 


Thus, the properties of S c(k ) follow from those of the Jost func¬ 
tion Si; here, we shall make use of the fact that if Tt(k) has a 
zero located at k - p (in the complex plane of the momentum), 
then the S -matrix has a pole at k - p and a zero at k = -p. 
Note that for any partial wave ( > 1, two poles will approach 



Temperature (K) 


Figure 6: Temperature dependence of rate coefficients for reaction (a) and elas¬ 
tic scattering (b) corresponding to the cross sections shown in Fig. [5] 

the threshold and collide at k — 0, while for s-wave ((' = 0) only 
a single pole at a time may reach the threshold. In this paper we 
analyze the case of s-wave scattering, for which it is convenient 
to factor out the effect of the pole at k — p and its accompanying 
zero at k - -p. Thus, we can write 

S(k) = ^S(k), (8) 

p — k 

and we assume the background contribution S ( k ) is a slowly 
varying function of k. Note that we shall omit the subscript t, 
as we only discuss the case of C — 0. 

As is well known 0], for purely elastic scattering in partial 
wave ( = 0, a pole p which is sufficiently close to k = 0 will be 
located on the imaginary axis; thus, we can write 

P = iK, ( 9 ) 

with k real valued; when k > 0 the pole corresponds to a bound 
state, while for k < 0 it corresponds to a virtual (or anti-bound) 
state. Since we are interested in near threshold resonances, we 
shall focus our attention on a short segment of the trajectory of 
the pole along the imaginary axis near k — 0. From Eqs. © and 
([8]) it is clear that the resonant contribution S Tes (k ) = jdi = j£±| 
is explicitly unitary for real values of k. The full S matrix is 
also unitary for real k, and thus the background part has to be 
unitary, and both of them can be written in the usual fashion, 
i.e., S = e 2 ’ 6 and S = e 2 ' 5 . The full phaseshift 6 is a sum of the 
background phaseshift 6 and the resonant contribution, 

6(k) = 6{k) - arctan(£) . 

Hence, the elastic cross section, cr elast = sin 2 5, reads 

elast// \ (k COS d-Af Sind) 2 

^ ^ k 2 k 2 +K 2 
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Figure 7: Mass dependence of the reaction rate coefficient for F + H2(0,0) —» 
HF + H at T = 1(T 6 K. 



Figure 8: Energy dependence of the reaction cross section for F + H2(0,0) for 
different masses of H. 


At low energy we use the approximation sinb « 5{k) ~ —ka 
(and cos 5 « 1) to further simplify the expression above: 


o- elast (/t) 


^ (1 +KO)- 2 

= 4ra j 


( 10 ) 


where a is the background scattering length, and a - a + - the 
full scattering length. We remark that the Lorentzian expres¬ 
sion in the equation above is a function of the momentum k. 
Thus, although Eq. ( fTTTb is similar to the Breit-Wigner formula, 
we recall that the latter has the familiar Lorentz-type energy 
dependence. Moreover, since the lineshape in Eq. CS is cen¬ 
tered on the threshold (k = 0) itself, only the k > 0 half can 
be probed in elastic scattering, which is specific to the case of a 
near threshold resonance for partial wave £ — 0 in the entrance 
channel. 

When the resonant pole crosses the threshold, the scattering 
length diverges: 


a-a + K 1 -> ±oo. (11) 

Hence, the behavior of the cross section will be affected dramat¬ 
ically; indeed, when k is vanishingly small, the denominator in 
Eq. m will give rise to two different types of behavior inside 
the low energy domain, depending on the relative magnitude of 
k and k. When k <sc |/c| we recover the familiar result for the 
Wigner regime, 


cr elast (k) 


k-> 0 


471a 1 , 


while for the remainder of the low energy domain, k becomes 
negligible when k » |/c|, and we obtain 


cr 


,elast 


(k) 


k»\K\ 


4 7t 


(1 + Ka) 1 

ie- 


4n 

¥’ 


which we refer to as NTR-regime behavior. From our discus¬ 
sion, it is clear that the extent of the NTR regime is specified 
approximately by the inequalities |/c| ;$ k IS \a\~ l . 

We emphasize that although the Wigner regime can be dis¬ 
placed towards the extreme low energy domain, Wigner’s thresh¬ 
old law is always recovered when k —» 0, except in the critical 


case of k = 0, which corresponds to a bona fide zero-energy 
resonance. Only in this special case, could the Wigner behav¬ 
ior be lost, and cr elast (k) would diverge at k - 0, since it would 
follow the NTR-regime behavior throughout the entire low-/: 
regime. As we shall see in the next section, the special case 
of a zero-energy resonance is only possible for purely elastic 
scattering, because any additional (inelastic or reactive) open 
channels—which are coupled, however feebly, to the entrance 
channel—will push the resonance pole away from the imagi¬ 
nary axis; hence, the pole will be kept away from k - 0, but 
still possibly close enough to produce a significant resonance 
effect. 

4.2. Pole factorization; many channel case 

We consider first the reaction cross section, followed by the 
elastic cross section. 

4.2.1. The reaction cross section 

We now extend the analysis of the single channel S matrix 
to the many channel case; namely, we focus on the diagonal el¬ 
ement S a(k) corresponding to the entrance channel. It is impor¬ 
tant to note that the factorization employed in the single channel 
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Figure 9: Energy dependence of the reaction rate constant for F + Fb(0,0) for 
different masses of H. 
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case in Eq. © also holds for the diagonal element correspond¬ 
ing to the entrance channel; indeed, the general properties of the 
S matrix ensure that a pole atk- p is always accompanied by 
a zero of S n(k) at k - —p |41J. Thus, for a NTR in the entrance 
channel /, we write 


Su(k) = ( 12 ) 

p-k 

The background contribution §u(k) = e 2,s ‘ is again assumed to 
be a slowly varying function of k. However, unlike the single¬ 
channel case, the background phase shift 6, is now a complex 
quantity (owing to the nonzero reactivity of the system). More¬ 
over, the pole p is no longer on the imaginary axis, as will be 
illustrated by our results, and can be explained by the fact that a 
quasi-bound state of the van der Waals complex near the thresh¬ 
old of the entrance channel can decay in all available (open) 
channels, and thus has a finite lifetime and a finite decay rate. 
Should the pole be located on the imaginary axis in the momen¬ 
tum plane, it would then be located on the real energy axis in 
the energy plane, implying a zero width (decay rate T = 0) and 
thus an infinite lifetime of complex, which would contradict the 
assumption of nonzero reactivity. 

Our results will show that Re(p) is always negative; thus, 
we write explicitly: 

p — —q + ik, 

with q > 0. The nonzero value of q stems from the nonzero 
reactivity (inelasticity) of the system, as argued above. Hence, 
the eigen-energy E p = p 2 /2p associated with the pole is now 
a complex quantity, with R e{E p ) = (q 2 - K 2 )/2p and Im(£,,) = 
-qx/p. The algebraic sign of k will determine if the pole repre¬ 
sents a virtual state or a quasi-bound state of the van der Waals 
complex; indeed, just as in the single channel case, k > 0 corre¬ 
sponds to the latter, while k < 0 to the former. We remark that 
the decay rate T = -2 Im(£’ p ) = 2 qK/p is positive for a quasi¬ 
bound van der Waals complex, while for a virtual state we have 
T < 0 (unphysical). 

At low energy, we employ the usual parametrization of the 
background contribution, <5/(k) » - ka where the background 
scattering length a, = a, - ifi, is now a complex quantity. Using 
similar parametrizations, the resonant part and the full S matrix 
element in Eq. (fl2l) can be linearized when k —» 0. Thus we 
have S™(k) = g * 1 - 2 ia™k, and S n (k) = S“(*)3„(*) « 
1 - 2 icijk, and we can extract the resonant contribution to the 
scattering length, aj es = j — * 2 ~ , and the full scattering length 
aj = aj + aj- es = a, + A Omitting the channel subscript, the real 
and imaginary parts of the scattering length read: 

Re(a) = a - a + (13) 

q - + k~ 

-lm(a)=/3=p+^-^. (14) 

q - + k~ 

We now return to the unapproximated expression of the 
resonant contribution in Eq. S res (k) = and the full 

S matrix element, S ( k ) = e ~' 2,ka , which are valid for the 

entire low-energy domain, namely k \d\. Making use of 



Figure 10: Low energy behavior of the total reaction probability for m = 1.042 u 
in H 2 +CI. The horizontal dashed line marks the unitarity limit. Inset: linear- 
scale plot of P react (k ) for different masses (as indicated). The data points from 
the full (exact) computation are shown as symbols, while the continuous lines 
were obtained by fitting using Eq. Gq). 


|5(k)| 2 \e 2,ka \ 2 = e 4k P, and substituting Eq. (ITH i in Eq. ©, 
the reaction cross section reads: 


£j_react 


7r 2 e- 2 ^ 
k 2 |p - k\ 2 


[(k 2 + |p| 2 ) sinh 2 kfi + 2gkcosh 2k/))j, 

(15) 


Assuming low reactivity, i.e.,/) —> 
simplified: 

cr rcact (k) * 47 t— 
k 


0, this result can be further 


p-k 


06) 


with [i given in Eq. m. Finally, if one neglects ft entirely, one 
can substitute /3 ~ /T es = q\p\~ 2 in the equation above, which 
can be recast as 


4-7T 

<x rcact (k) » — 


q 


k (q + k) 2 + k 


■2 ' 


(17) 


This expression only contains two parameters (q and k) which 
are sufficient to describe a prominent NTR (when is small). 
Note that the full expression given in Eq. lfl5l > includes the back¬ 
ground contribution ft as a third parameter, which allows for 
slight deviations from the simple k-dependence in Eq. dlTb . 

Recall that we are interested in the case when \p\ is very 
small, such that the energy associated with the pole, \E p \ = 
4j}h 2 \p\ 2 , is well within the ultracold domain. Consequently, 
there will be two different types of behavior at low k, depend¬ 
ing on the relative values of k and \p\. Indeed, for k <s \p\, we 
recover Wigner’s threshold law. 


^eact -, An P_ 

k-> o k 


( 18 ) 
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while for k » \p\, Eq. (fTTb reduces to 


tT™* - > 4 4 

k »|p| k 3 


(19) 


which we named NTR-regime behavior. 

We emphasize that the simple expression in Eq. ( fTTb cap¬ 
tures both types of behavior, including the transition between 
the two regimes, and we can employ it to fit our computed 
data. Specifically, we use a more convenient quantity than the 
cross section shown in Fig. 0 a), namely the reaction probabil¬ 
ity, H react = 1 - I.S’,,1 2 = § cr react (k), which is bound between zero 
and one. Hence, our fitting procedure employs the expression 


\k) = 4 q 


C q + ky 


( 20 ) 


which is equivalent with Eq. G3- The asymmetrical profile of 
P react (k) is clearly illustrated by the linear-scale plot shown as an 
inset in Fig.[T0l which contains results for various masses, rang¬ 
ing from the true mass of hydrogen ( m — mu) to m = 1.042 u. 
The results for one of the strongly resonant cases (m = 1.042u) 
are also shown on a log-log scale (see the main plot in Fig.lTOli. 
which emphasizes the good agreement between the numerical 
results and the simple analytical expression in Eq. (l20l) . Also, 
the two power-laws in Eqs. (ITSb and (fl~9b are easily recognized 
on the log-log plot in Fig.|T0j indeed, we have P K ' dCl (k) ~ k in 
the Wigner regime, and P redct (k) ~ k~ l in the NTR regime. 

From Eq. (l20l > it follows that P react (k) reaches its maximum 
value at k max = \p\, which marks the transition between the two 
regimes; in terms of energy, the maximum of f >react is located at 
Emax = \E p \ = T^iq 2 + k 2 ). Although it is tempting to interpret 
the location of the maximum of P react (A;) as the position of the 
resonance, one has to resist this temptation. Indeed, based on 
Eq. (l20b . we suggest that it is more appropriate to assign k - -q 
as the position of the resonance, instead of k - k max = \p\. 

We remark that the unitarity limit (_P react = 1) can only be 
reached in the special case with k - 0, when the pole is situated 
on the negative real axis at p — -q, and the diagonal element 
of S corresponding to the entrance channel has a zero on the 
positive real axis at k - -p - q. Namely, we have S ,,(q) = 0, 
see Eq. (IT2b . which guarantees that P redct (q) = 1. Figure Qj] 
shows the pole trajectory in the complex plane as a function 
of the parameter m, and makes it apparent that the trajectory 
p{m) crosses the real axis when m reaches a critical (resonant) 
value m - m res . We note that in the special case of m = m res , 
when the unitarity limit P ledct (k) = 1 is reached at k - q(m les ), 
the NTR regime extends to the lowest possible energies into 
the ultracold regime: k ;$ \p{m K% )\ — q(m K% ). When m moves 
gradually away from m les the NTR regime shrinks gradually, 
and it spans a narrower range of energies (as seen in Figs. 5, 
6, 8, and 9). Eventually, the NTR regime disappears when the 
pole p{m) is no longer in the vicinity of k — 0. 

The pole parameters q — -Re(p) and k - \m(p) are treated 
as fitting variables, and they are thus obtained from the com¬ 
puted results. Note that the value of /S can only be extracted 
from the fitting procedure if the full expression in Eq. O is 
employed instead of the simpler Eq. (EJ; we also remark that 


Im[k] (a.u.) 



Figure 11: Pole trajectory in the complex plane of the momentum for H 2 + Cl. 
The inset shows the explicit mass dependence of q = -Re[/?] and k = Im[p]. 
Note that q(m) is nearly constant, while K(m) has an almost linear behavior. 
Also note that |/c| increases considerably on either side of the resonance, such 
that q |/c| (except near the very center of the resonance, where k « 0). 


when p is negligibly small compared to the resonant contribu¬ 
tion /l res = q(q 2 + icy 1 , see Eq. ( ITTb . then it can no longer be 
extracted reliably from the fitted data. 

We emphasize that we conducted our analysis in the com¬ 
plex plane of the momentum k, rather than the complex energy 
plane; to appreciate the unsuitability of the latter, we simply 
point out that the complex energies associated with both the 
pole (k = p) and the zero (k = -p) are identical. 


Er = 


n 2 (±p ) 2 

2p 



( 21 ) 


with s p = -r(q 2 - K 2 ) and E = —lik. We recall that one must 
have q > 0 for a physical cross section, see Eq. (fl7b . and thus 
we have sgn(T) = sgn(/c). Hence, for a quasi-bound vdW- 
complex (k > 0) we have T > 0, while for the case of a virtual 
state (k < 0) we have T < 0 (unphysical). Finally, we remark 
that for an ,v-wave NTR, scattering alone cannot distinguish be¬ 
tween shallow bound and anti-bound (virtual) states, because 
cr rcacl (/:) does not depend on sgn(A-), see Eq. ©. 


4.2.2. The elastic cross section 

For the sake of completeness, we also give the expression of 
the elastic cross section. Employing again the parametrization 
§ = e ~ 2lak e~ 2 P k in Eq. (fl2l) . which is then substituted in Eq. ©, 
we obtain 

Ti 2e~ 2k P 

^eiast _ -_ |(£ 2 + |p| 2 )cosh 2kfi + 2kKsin 2ka 

k~ |p - k\ z L 

+ ik 2 - |/?| 2 ) cos 2 ka + 2kq sinh 2k/)J. 

Under the assumption of low reactivity, p —> 0, which is a re¬ 
quirement for having pronounced NTRs, we obtain 


cr 


elast 


4zr|a| 2 


P 

p — k 


2 


where a — a - ip is the full scattering length. Although the ex¬ 
pressions for the elastic cross section given above contain four 
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parameters (a, /3, q and k), only two of them are needed when 
the resonance pole is very close to the threshold (p ~ 0). In¬ 
deed, if the pole contribution is dominant, we can neglect a, and 
use a « a res — j } in the equation above, which we rewrite as 


cr 


elast 


4 71 

(,q + k ) 2 + k 2 


( 22 ) 


Note the similarity between this result for the elastic cross sec¬ 
tion and Eq. (ITTI i for the reactive cross section; thus, in the limit 
k —> 0 we recover again the familiar Wigner behavior. 


cr 


elast 


-> 47r| a\ 2 , 

k-> 0 


and for the NTR regime we obtain 


cr 


elast 


-> 

k^>\p\ 


4 7T 

¥' 


5. Jost-function approach 

We now describe an alternative approach, based on the use 
of the Jost function, to obtain the anomalous k ~ 3 scaling of the 
reaction cross section for the NTR regime described in the pre¬ 
vious section. We shall first give the key concepts for the case 
of single channel and short-range interaction V(R), and gener¬ 
alize the results to the multi-channel case afterward. 


5.1. Single channel case 

The regular solution of the radial equation, 

d 2 

—<P k = \2liV(R) - k 2 ]<p k (R), 

is uniquely specified by initial-value type conditions: (pk - 0 
and -jx<pk = 1 at R = 0. Recall that we are discussing only the 
case of s-wave (l = 0) scattering; for higher partial waves, the 
initial boundary condition reads: <pk,e(R ) ~ R e+l . The asymp¬ 
totic behavior of the regular solution can be written as 

<l>k(R) — - * — \A(k) sin kR + B(k) cos kR], (23) 

k 

If k is real valued, and assuming V(R) is real, then we can 
choose tpk to be real; hence, A(k) and B(k) are real valued. Re¬ 
casting Eq. (1231) in terms of the free solutions exp (±ikR), we 
have 


MR) Yk t (A _ lB)e ~ m ~ {A + lB)e+,kK \> ( 24 ) 

and following the convention used in Refs. 00, we define 
the Jost function: 


F(k) = A(k) - iB{k). (25) 

Conversely, if one introduced F in the usual manner, i.e., by 
starting with the asymptotic behavior of (pk written as 

M R ) Yk \ T(k ) e ~ lkR - r*(k)e + ‘ kK \. (26) 


one could identify A = RefF) and B — -Im((F), provided that 
both k and ViR) are real valued. Note that if one allows for 
a complex valued k, see 0, a fully general definition of the 
Jost function can be obtained by replacing F*(k) with 'F(-k) 
in Eq. (l26t . In this article we use real and positive k. unless 
otherwise specified; when needed, we can allow k to be com¬ 
plex without any difficulty, because we restrict ourselves to the 
case of short-range interactions. We remark that our analysis 
of the S matrix in the previous section (done in the complex 
k plane) would be difficult to justify rigorously for the case of 
long-range interactions. However, as we shall see in this sec¬ 
tion, the approach based on the Jost function (for real argument 
k) circumvents this obstacle, because it makes it possible to ac¬ 
count for the NTR phenomenon without referring to the poles 
of the S matrix in the complex plane. Moreover, when it is pos¬ 
sible to define the Jost function for complex k, it is easily seen 
that the two approaches are equivalent. 

Regarding our preference for analyzing A(k) and B(k) in¬ 
stead of F(k) itself, we point out that although the two solutions 
xf(R) ~ e ±,kR are linearly independent in the strict sense, they 
will become linearly dependent in the limit k —> 0 (from the 
numerical point of view) simply because e~ lkR ~ e ,kR ~ 1 when 
kR <s 1. Thus, it is desirable to replace xf with a set of solu¬ 
tions which are more suitable at low k, namely fi : (R) ~ sin (kR) 
and gk(R) ~ cos (kR). Since we assume that M fk, gk and x\ are 
exact solutions of the radial equation, the matching conditions 
(1231) and ( 1261 ) can be recast as equalities: kipk = Afk + Bgk, 
and 2ik<pk = F*xt ~ Fx l; respectively. As is well known, 
the Jost function can be expressed as F — W\xk, <PkL where 
W[u, v ] = uv' - u'v denotes the Wronskian. Similarly, for A and 
B we have A = W[gk,(pk\ and B — -W[fk,(pk], which can be 
used in practice to compute A(k) and B(k). 

We now shift our attention to the (normalized) physical so¬ 
lution \pk{R ), which is regular at R = 0. Hence, the physical 
and regular solutions only differ by a multiplicative constant, 
which can be easily found by comparing the as ymp totic bound¬ 
ary condition obeyed by the physical solution |43 m. 

MR) > e iS sin (kR + 6) = (e~ ikR - e 2is e ikR ) , (27) 

to the asymptotic behavior in Eq. (l26l >. Indeed, we find 

MR) = -MR), (28) 

F{k) 

and we can also identify the S matrix, S — e 2 "' - F*F~ l ■ We 
recall the well known relationship between the phase shift 6(k) 
in Eq. (1771) and the Jost function; namely, S = — arg (F), which 
can also be written as tan 6 - B/A, see Eqs. (I23ll25l[27l >. Also 
note that if we make use of F — \F\e~' 6 , then Eq. (l25l) can be 
expressed as A — \F\ cos 6 and B = \F\ sin 8. 

We emphasize that scattering resonances will affect strongly 
the k dependence of the normalization constant k[F(k)] in 
Eq. (l28l >. Special attention has to be paid to the denominator 
F — A — iB, which could vanish near k - 0, and thus produce 
a near threshold resonance. Hence, we proceed to analyze the 
low k behavior of the Jost function. Under the assumption of a 











short-range potential, it is possible to expand A(k) and B(k) as 
power series. Indeed, given that the regular solution, as well as 
cos (kR), are analytic in k 2 [440, and that sin (kR) is analytic in k 
and is an odd function of k, we can write: 


Thus, the elastic cross section reads 


cr elast (£) 


4tt B\k) 

T 2 A 2 (k) + B 2 (k) 


( 31 ) 


A(k) = A 0 +A 2 k 2 +A 4 k 4 + ■■■ (29) 

B(k) = Bik + B 3 k 3 + B 5 k 5 + ■■■ (30) 


For a finite-range potential, these power series can be rigorously 
deduced in a straightforward manner, and their validity can also 
be extended to the case of short-range potentials. However, for 
long-range potentials, e.g., V(R) - -C n /R n for R —> oo, A(k) 
and B(k ) are no longer analytic functions; their k dependence 
is rather complicated, as shown by Willner and Gianturco [451. 
Nevertheless, even in the case of long-range interactions, both 
A(k) and B(k), and hence T(k) are generally well behaved, al¬ 
beit nonanalytic; thus, we subscribe to the point of view advo¬ 
cated by Willner and Gianturco [45f], that the low-/: expansion 
for the Jost function is simpler, more elegant, and more practi¬ 
cal than the corresponding expression for tan 6(k) or the famil¬ 
iar effective range expansion for k cot 6(k). We argue that one 
should employ the low-/: expansions for A(k) and B(k) of the 
Jost function Tik) — Aik) - iB(k) rather than low k power series 
expansion of the phase shift (and all other physical quantities, 
such as the cross section). For example. 


tan<5(k) = 


B(k) 

A(k) 


B\k + B 2 k 3 + • • • 
A 0 + A 2 k 2 + ■■■ 


should be kept as a fraction, rather than attempting to expand 
it consistently to 0(k 2 ) or higher. Indeed, such an expansion 
would become unsuitable if Ao were vanishingly small, as its 
validity would be restricted to an extremely narrow range, namely 
k 2 <s IA 0 /A 2 I, while the expression above remains valid and 
useful for a much wider range. Alternatively, we have 


k cot d = 


Ao + A 2 k~ + ■ ■ • 
Bi + B 2 k 2 + • • • ’ 


which could be expanded as a power series to give the familiar 
effective range expansion. 


1 4 . c ^ ^eff ,7 

k cot <5 w- 1 - lc, 

a 2 


with a = -Bi/Aq the scattering length, and r e g = 2(A 2 /B\ - 
AqB 2 /B 2 ) = (2/Bi)(A 2 + B 2 /a ) the effective range. When Ao 
vanishes, we have a —> ±00 and r s g —» 2A 2 jB\, and the effec¬ 
tive range expansion remains valid, taking a simplified form: 
k cot <5 = \r e gk 2 . However, if B\ is vanishingly small, we have 
a * 0 and r e f[- —> + 00 , and thus the effective range expansion is 
rendered useless. 

We now analyze the elastic cross section in terms of the Jost 
function: 

cr eUs \k) = ^\l-Sm\ 

where the S matrix is expressed as 


In the limit k —* 0, and provided that A(0) A 0, we recover the 
well known Wigner law for s-wave elastic scattering, which is 
simply the lowest order approximation. 


cr 


elast 


-» 4-7T 

/ t->0 



= 47Tfl 2 . 


The domain of validity for Wigner’s threshold law defines the 
so called Wigner regime; its extent is governed by the lowest or¬ 
der coefficients Ao and B\ in the Jost function expansions (l29l 
m. When neither Ao nor B\ is small, the Wigner regime will 
cover the entire low energy domain. Conversely, if either Ao 
or B\ is vanishingly small, Wigner’s threshold law will be re¬ 
stricted to a narrow domain near k = 0, because the higher or¬ 
der terms in ( l29l l30l > become dominant as k increases, and the 
^-dependence of the cross section in Eq. <EQt changes dramati¬ 
cally. 

Here, we are especially interested in the case of vanish¬ 
ingly small Ao, which corresponds to a NTR, and we estab¬ 
lish the domain of validity for Wigner’s law as follows: in 
Eq. (ED i s necessary that both A(k) and B(k) be replaced 
by Eqs. ( l29l l30l > truncated to the lowest order, i.e., we require 
k «c IA 0 A 2 1 1 1/2 and k «c 1 1 1/2 . The former inequality is 
stronger than the latter, since B\ is typically finite when Ao —> 0. 
Moreover, the B 2 (k) term in the denominator of Eq. ( OH has 
to remain negligible; hence, we impose yet another inequality, 
k <s IAoBj 2 ! = |a| _1 , which is much stronger than either of 
the two inequalities above, and thus dictates the extent of the 
Wigner regime. 

Since the Wigner regime is displaced towards extremely 
low energies when Ao ~ 0, we need to analyze the behavior 
of the elastic cross section for the remainder of the low energy 
domain. Thus, we consider k outside the Wigner regime, but 
still inside the low energy domain, such that Eq. POl) can be 
truncated to the first term, and Eq. ( l29t to the second term. As 
k increases gradually, Ao becomes negligible compared to both 
A 2 k 2 and B\k. Moreover, only the latter term survives, as we 
have \B\\k » IA 2 IZ : 2 throughout the low energy domain. Con¬ 
sequently, the term A 2 (k) in the denominator of Eq. (OH can be 
neglected entirely, and we recover the NTR-regime behavior, 

*-(*)*£■ 


In summary, using the Jost function description, the elastic 
cross section can parametrized in a simple way at low energy; 
namely, according to our discussion above, we keep only the 
relevant terms in Eqs. (l29l[30l > to approximate Eq. OH as 


cr 


,elast 


( k ) = An 


B\ 


A 2 + B\k 2 


= 47T 


1 


k 2 + (1/a) 2 ’ 


(32) 


S(k) = e 2is 


A + iB 
A-iB' 


where a = -B\ /A (l is the full scattering length. The simple 
expression above makes it clear that the the gradual transition 
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from Wigner’s threshold law to the NTR-regime behavior takes 
place near k ~ |a| -1 , where the two terms in the denominator 
are comparable. We remark that Eq. ( l32l > is in agreement with 
Eq. (fTOl) obtained in the previous section using the pole factor¬ 
ization of the S matrix, and we shall discuss their equivalence 
in Sec. [6] 

5.2. Multi-channel case 

In the case of N coupled channels, the regular solution 
is an N x N matrix with elements satisfying the system of 
coupled radial equations, 

-^<D(k, R) + K 2 0(k, R) = ^V(R)®(k, R), (33) 

with k s (k\, ki, • ■ ■, k;v) an d k n = y/2/u(E - s„)/H 2 the mo¬ 
menta corresponding to channel threshold energies s„. For sim¬ 
plicity we consider all channels open, with all k„ real valued 
and k n >0. K is a diagonal matrix with elements k„6 nm , while 
V(R) is a full matrix containing the couplings V„ m (R). Note that 
the diagonal elements of V also include centrifugal terms of the 
formf’„(f’„ + 1)R 2 . 

The regular solution is uniquely specified by initial-value 
type conditions at R — 0, 

it—>o , 

Kn(R) -* 6 nm R e " ■ 

Assuming that the coupling matrix V is real, we extract two real 
valued matrices (A and B ) from the asymptotic behavior of the 
regular solution, 

ch -> K"' _1 (jA + nB), 

where K = diagfk,/" -1 ), and the diagonal matrices j and 
n contain the free solutions, i.e., the Riccati-Bessel functions 
j((k„R) and ii((k n R) for each channel. Note that, although the 
asymptotic behavior in the last equation was written in matrix 
notation, the coefficients A nm and B nm (i.e., the matrix elements 
of A and B) are extracted independently for each pair of indices 
(n, m). Indeed, for each component (n) of any column solution 
(m) contained in <X>, we write the asymptotic behavior explicitly, 

p , i /?—» co 

kn * A.nmjt n (fcnR) B nm Yli n ih n K) • 

Employing this matching condition for <f> nm (together with a 
similar equation for dcf> nm /dR ) at R - R max in the asymptotic 
region, we obtain 

A nm = ~k e n a W[tf> nm , n { J (34) 

Bnm =+k e n n W[(p nm , j( n \ (35) 

where W[... ] stands for the Wronskian. The Jost matrix is de¬ 
fined in terms of the matrices A and B as 

F = A-/B. (36) 

which is a direct generalization of Eq. (l25l >. The factor k { " +x 
in the equations above was used for convenience, such that the 


Jost matrix is well behaved when k n —* 0. In this paper we 
restrict our discussion to s-wave scattering, i.e., i — 0 in the 
entrance channel; note that only the centrifugal term in the en¬ 
trance channel is relevant for analyzing the threshold behavior. 
The values of £„ in all other channels can be arbitrary (as al¬ 
lowed by the specifics of any given scattering problem). How¬ 
ever, we set £„ = 0 in all channels to simplify our notation; thus, 
the asymptotic behavior of the regular solution reads 

/?—» oo 

k-ntynmiR) * sin (k n K) + B nm cos (k n R), 

which is identical to the single channel version, see Eq. dm 
Equivalently, the asymptotic behavior of O can be written in 
terms of Jost matrix elements F nm = A nm - iB mn as 

—kd> (R) R > e~ ik " R F - e +ik '< R p’* 

. 'WzY'nmv'V 7 c r nm " r nm’> 

where F* m = A nm + iB„ m is the complex conjugate of F„ m . 

The normalized (physical) solution V can be written in terms 
of the regular solution <[>, 

^(R) = <D(R)F -1 K, (37) 

and it has the well known asymptotic behavior, 

^nm(R) ^ \e- ,KR 6 nm - e +lk ” K (k„- |/2 5„ m k;„ /2 )|, (38) 

where S are the elements of the S matrix, 

s = kr 1/2 (F)*(Fr I K 1/2 . (39) 

Equations 071 ) and ( l39l > contain the inverse of the Jost matrix, 
which we write explicitely as 

F 1 = [Co f ( p )] r , (40) 

det(F) 

with [Cof(F)] r the transpose of the matrix of cofactors of F. 
and det(F) the determinant of F. 

As in the single-channel case, we focus our attention on 
the denominator, i.e., the determinant of the Jost matrix. We 
are interested in the situation when det(F) is vanishingly small, 
such that the scattering cross sections are resonantly enhanced; 
if T„ m are the elements of the T matrix, T = 1 - S, we have 

cr(n <- m) ~ -^-\T nm \ 2 ~ \ det(F)|“ 2 . 
km 

Although the determinant of the Jost matrix cannot vanish on 
the real axis |[46i], it can almost reach zero. Also, recall that we 
only consider the case of a scattering problem without closed 
channels, thus eliminating the possibility of Feshbach resonances 
Therefore, by extrapolation from the single-channel case, the 
only remaining possibility is that of a potential resonance at low 
energy; such a resonance would correspond to a quasi-bound 
state near the threshold of the channel with the highest energy 
asymptote, which we take as our entrance channel (n = 1). For 
clarity, we assume E > s\ > £2 > ■ • • > e n > ■ ■ ■ > sn- 
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When det(F) * 0 at a channel threshold, the factor | det(F)| 2 
in the equation above will be responsible for a resonance en¬ 
hancement, and will also affect the threshold behavior of cross 
sections. We now give a detailed analysis of a NTR in the en¬ 
trance channel, n — 1, and we start by using a cofactor expan¬ 
sion of the determinant along the first row, 

N 

det(F) = ^ F lm Ci„„ (41) 

m= 1 

where C\ m is the cofactor corresponding to the element (1, m). 
Next, we extract Cn outside the sum, and we isolate the first 
term (i.e., the diagonal element Fn) to recast det(F) as 

det(F) = C n (F 11 +/ 11 ), (42) 

where 

fu=-^yF lm C lm . (43) 

Cl1 

Separating the real and imaginary parts of the last factor in 
Eq. (l42l >. namely Fn + fn — JK - iS, and defining 

D = \F n +fu\ 2 = AR 2 +S 2 , 

we can write 

I det(F)| 2 = \C n \ 2 \M-i!B\ 2 = \C n \ 2 D. (44) 

Under the assumption that fn is suitably small, and pro¬ 
vided that Cii # 0, only the factor D ^ \Fn\ 2 - A 2 U + B 2 U can 
be responsible for NTRs. In order to clarify these assumptions, 
we recall that, as we explained in Sec. [4] due to the coupling of 
the entrance channel with other open channels, the resonance 
pole will be pushed away from the threshold; thus, in the limit 
of strong coupling, the existence of a near threshold resonance 
will become highly unlikely. Conversely, the appearance of a 
prominent NTR requires the entrance channel be weakly cou¬ 
pled to all other open channels. To simplify, one can use a 
coupling strength parameter A to redefine these particular cou¬ 
plings (Vi„ = V n i —> AV n i) which will become vanishingly 
small when A —» 0. Hence, we assume the regular solution 
matrix <J) is almost block diagonal, with a small lxl block 
containing (f>n and a large (N — 1 ) x(N - 1) block correspond¬ 
ing to all other open channels. The remaining elements, i.e., 
with m + 1 in the first row, and <p n \ with n + 1 in the first 
column are vanishingly small when A —> 0. Consequently, from 
Eqs. (l34l1351 l36l > we have (for n + 1) F n \ = <9(T) when A —> 0, 
and thus C\ m = 0(A) for m + 1. Note that Cn is independent 
of /I in zeroth order, and hence we assume Cn - 0( 1). Finally, 
we also have F\ m = 0( A) for m + 1, and from Eq. (l43l > we have 
fn = 0(A 2 ) when A —> 0. 

Regarding the caveat of an accidental vanishing of Cn near 
the threshold of the entrance channel, recall that in the limit 
A —> 0 the channel n — 1 becomes decoupled from the remain¬ 
ing N — 1 channels, and we have lim,i_>oCii = det[F (W_ •*], 
where F (,v_1) is the Jost matrix for the scattering problem with 
N — 1 channels. Thus, det | F' 11 ] cannot be vanishingly small 
for collision energies E > s\, since E is very high above £2 


(which is the highest threshold for the scattering problem with 
N — 1 channels). Indeed, as we assume e\ - £2 is much larger 
than the generic low-energy scale of the scattering problem; 
thus, it is safe to factor out Cn = <9(1) as a background contri¬ 
bution, as we already anticipated in the equations above. 

Denoting k - k\, we now focus on the k dependence of D(k) 
in Eq. (l44l >. which stems from the behavior of Fn(k) and fn(k). 
We emphasize that all quantities are analyzed as functions of 
the single variable k. All other momenta, i.e., k„ for channels 
n = 2,3,... ,N, can be easily expressed in terms of the entrance 
channel momentum: using k 2 = 2 fi(E - s n ), we write 

k n — Va„ + k 2 , 

with A„ = 2/r(si - £„). In the ultracold limit (k —> 0) we have 
k 2 «c A„ for all n + 1, and we use a power series expansion: 

k n (k 2 )* VA ~ + 0(^). 

Moreover, the regular solution <X> is a function of k 2 (rather than 
k itself), and we deduce from Eqs. (l34l [35l > that except for the 
first row, all matrix elements A nm and B nm are also functions of 
k 2 . Thus, when k —> 0 we have 

F„ m (k 2 ) = F nm ( 0) + 0(k 2 ), for n+\. 

Similarly, the cofactors corresponding to the first row of F can 
be written as 

C hn (k 2 ) = Ci„,(0) + 0(k 2 ), for all m, 

which ensures that the ratios c m = Ci m /Cn in Eq. (1431) can be 
expanded in a similar fashion: 

Cm(k 2 ) - Cm(0) + 0(k 2 ). 

Rewriting Eq. (S3 as 

fu = ^ F lm c m, (45) 

m +1 

and separating its real and imaginary parts, we have 

Re(/n) = X \AUmccJk) + B hn (k)lmc m (k)], (46) 

m±\ 


Im(/n) = ^ \A Un (k)\mc m (k) - B lm (k)Rec m (k)\. (47) 

nvk 1 

In these equations, special attention has to be paid to A\ m (k) and 
B\ m (k). From Eqs. (I34ll35l ) we obtain the expansions: 

A hn (k) = A,„M + k 2 A {2) m + 0(k\ (48) 

BUk) = kB^+^Bfl + Oik 5 ). (49) 

Since the series of B\ m (k) contains odd powers, and c m (k) in 
Eq. S3 af e complex quantities, the power series in Eqs. (l46l 
1471) contain both odd and even powers of k, 

Re(/„) = a(0) + ka 0) + k 2 a (2) + 0(k 3 ), (50) 

Im(/n) = b(0) + kb (1) + lrb (2) + 0(k 3 ). (51) 
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Thus, the power series of m(k) and B(k) will also contain all 
possible terms, i.e., both odd and even powers: 

(k) = Jlo + km i + k 2 m 2 + ■ • • , (52) 

B(k) = ®o + kS i + k 2 B 2 + ■■■ (53) 

The coefficients m v and B v will be expressed in terms of 4‘’' 2 , 
a (v) and b (v \ Indeed, recall that m = Re(Fn + fn) = 
A n + Re(/n) and B = -Im(F n + f n ) = B n ~ Im(/n). From 
the equations above, we obtain for the coefficients up to order 
v = 2, 


mo = mm 
m x 
m 2 

and 

So = S(0) 

51 

5 2 


a(0)+A n (0) 




-b( 0) 

-/> (1) + B { " 

-b (2 \ 


Under the assumption of weak coupling (A —> 0), we have 
/n = 0(A 2 ) —* 0. Thus, all and /? (v) are negligibly small; 
however, their presence in the equations above is necessary, be¬ 
cause they limit the effects caused by near threshold resonances. 
Indeed, note that mo and Bo cannot vanish simultaneously; in 
other words, one cannot have: An(0) = 0 and a( 0) = 0 and 
b( 0) = 0. Were that the case, IXk) and hence det(F) would 
vanish at k — 0, which is impossible (as we mentioned earlier). 
Nevertheless, when A is small, so are a{ 0) and /;(()). Also, as 
in the single-channel case, it is possible that An(0) « 0. Thus, 
although the quantity D(k - 0) = ml + B 2 cannot vanish, it 
could reach very small values (which is the signature of NTR). 

Finally, we analyze the threshold behavior of cross sections 
in the presence of a NTR; from Eqs. (I39ll40l > we obtain 


Recall D(k ) = m 2 (k) + B 2 (k), which is expanded at low k, 

D(k) = Do + kD\ + k 2 D 2 + ■■■ , (56) 


with 

D 0 = D(0) = [/7(0)] 2 + [a(0) + A„(0)] 2 , 

D\ = 2[-b(Q)B\ + a w m Q ], 

D 2 = B 2 + [a (1) ] 2 + 2[m 0 m 2 + b(0)b (2) ]. 

As argued above. Do cannot vanish identically; however. Do 
is small in the case of a NTR, and hence it is the most impor¬ 
tant parameter characterizing the NTR. Note that the first order 
term in Eq. (l56l > can be neglected when both the weak coupling 
(A —» 0) assumption and the NTR assumption (4 1 ] (0) ~ 0) are 
valid; indeed, we have D\ ~ —2b(0)B l2 ] > — O(A). Thus, even 
for a pronounced NTR, i.e., when Do approaches its minimum 
value. Do ~ D m m ~ [/?(0f] 2 = 0(A 2 ), the linear term will be¬ 
come negligible compared to Do for k <s Do/\D\ \ ~ |/7(0)/Z? < 1 1 1 ) |. 
In the ultracold limit (k —> 0) the quadratic term in Eq. (l56l) 
can be neglected as well; specifically, the requirement is k <s 
\JDo / D 2 . Making use of the approximations Do ~ [/A0j] 2 and 
Z >2 ~ [B ( 1 1 1 ) ] 2 , we have \JDo/D 2 ~ |fe(0)/B l 1 1 I , |. Hence, the strong 
inequality k <s • \JDo/D 2 defines the Wigner regime, since we 
can ignore both the linear and quadratic terms in Eq. (l56t . In¬ 
deed, we have D(k) ~ Do, and we recover Wigner’s threshold 
laws: cr 16301 ~ /r 1 and cr elast ss constant. 

The linear term in Eq. (l56t can also be neglected outside 
the Wigner regime, i.e., for the remainder of the low k domain; 
indeed, when k » |Z?(0)/fi < 1 1 1 ) |, we have k 2 D 2 » kD\ and also 
k 2 D 2 » Do, which ensures that D(k) ~ k 2 D 2 . Thus, discarding 
the linear term altogether in Eq. (l56l t. we simplify Eqs. ([541 and 
(l55l >. and we obtain 

^O.l^elast, react ^ - (5 7 ) 

Do + k 2 D 2 


S run — b n 


k m 2/ 2 j B n jC mj 


k„ det(F) 


We only consider scattering from the ultracold channel m = 1 
to all final channels n = 1,2,3,..., A; thus, using the notation 
k = k\, the cross sections read 


where the prefactors k° = 1 and k correspond to elastic and 
reactive scattering, respectively. We emphasize that this ex¬ 
pression is a good approximation throughout the entire low-/: 
regime, and despite its simplicity, it accounts for the combined 
Wigner and NTR-regime behavior. For reactive scattering we 
have 


cr(n <- 1) ~ ^|7’„ 1 (Z:)| 2 ~ 


I Tij B n jCij\~ 

kk n \Cu\ 2 D(k) 


£j_react 


k 1 , Wigner regime for k <s \jDo!D 2 
k~ 2 , NTR regime for k » \jDo!D 2 , 


Except for D(k), all other quantities are finite when k —* 0, and 
can be replaced by their values at k = 0, i.e., B n j( 0), Cij(0), 
Cn(0) and k n ( 0), which amount to an overall (n dependent) 
constant. This allows us to simplify the expression above; for 
reactive (or inelastic) scattering, n + 1, we obtain 


cr react (n <- 1) ~ 


const(n) 

D(k) 


k~\ 


(54) 


while for elastic scattering, we have k„ = k and B\j(k) ~ k, and 
the cross section reads 


(T 


.elast 


(i 


const' 


(55) 


while for elastic scattering, 


(T 


elast 


k°, Wigner regime 

k~ 2 , NTR regime, 


According to Eq. (l57l > there is a smooth transition between the 
two types of behavior (Wigner and NTR) which takes place 
around k = \JDo / D 2 , where all three terms in Eq. < 1 56 b are 
comparable. Thus, the linear term in Eq. (l56l) only plays a mi¬ 
nor role near the transition, while becoming negligible for both 
Wigner and NTR regimes. 

We remark that when the coupling strength A increases, and 
the system becomes more reactive, a (v) and b (v) will increase; 
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Figure 12: Comparison of reaction probabilities for the resonances of H 2 + Cl. 
The imaginary part of the scattering length (ft, see Fig. [3) is shown in the inset 
to help identify the resonances. The three solid curves correspond to the three 
critical values of m (for the first, second and third resonance, which are shown 
in red, green and blue, respectively). The dashed curves correspond to values 
of m near the first and third resonance. The arrows indicate the increase of 
reactivity, and the increase of ft and q, as m decreases. 


hence, Dq and D\ will become large (even if Ao ~ 0 still holds) 
and the NTR regime will be pushed to higher values of k, where 
higher-order terms will become dominant. Eventually, the sig¬ 
nature of the NTR will disappear, as suggested in Fig. [T2] which 
shows that for the lowest reactivity ft corresponding to the third 
resonance, the pole is located at very low k, and the NTR regime 
extends to a sizable range for larger k. The extent of the NTR 
regime is smaller for the second resonance corresponding to a 
larger /?, and smaller still for the first resonance with the largest 
reactivity /i. Thus, in the limit of strong coupling, the very ex¬ 
istence of a quasi-bound state of the atom-dimer complex be¬ 
comes unlikely, since such a complex will be short lived; conse¬ 
quently, NTRs cannot exist (or are less pronounced) for highly 
reactive systems. 


The background contribution F(k) is also analytic, and assum¬ 
ing it to be slowly varying within a small domain containing 
both k = p and k - 0, we can truncate its power series in the 
low energy domain: 


F(k) ~ F 0 + F\k. 


Note that for F(k) we employ a power series expansion around 
k = 0, which is more convenient. Hence, the S matrix reads 

T(-k) p + k F(-k) 

T(k) p-k F{k) ’ 

and we recognize the background contribution, see Eq. 0, 


S(k) = 


F(~k ) 

m ’ 


which is slowly varying near k — 0, due to the smooth behavior 
of F{k). Truncating the power series of S to first order, we have 
S(k) as 1 - link. The unitarity of § on the real axis ensures 
that the coefficient of the first order term in this expansion is 
purely imaginary, and we have used the customary notation, 
namely the background scattering length a. Thus, in the pole 
factorization approach, the full S matrix at low k reads 

S(k) * Ejtl(i-2iak). 
p-k 

We now compare this result with the equivalent expression 
obtained in the explicit Jost function approach, namely in terms 
of F(k) = A(k ) - iB(k). We emphasize that A(k) and B(k) are 
expanded in power series around k = 0 only, and we remark 
that for a short range potential, the power series of A(k) and 
B(k) remain valid when k takes complex values; moreover, the 
expansion coefficients A„ and B„ in Eqs. (l29l [30l > remain real 
numbers. Therefore, if a zero of the Jost function (i.e., a pole 
of the S matrix) is located on the imaginary axis, p = ik, then k 
is a real valued solution of the equation 


6. Equivalence between different approaches 


Aq + B\k - A 2 K 2 - B^k 2 + ■ • • = 0, 


In the case of purely elastic scattering (single channel), the 
equivalence of the two approaches is based on the simple rela¬ 
tionship between the S matrix and the Jost function, see Eq. (0. 
The expressions 0 and ( I29ll30l) merely give different parametriza- 
tions for the S matrix, and hence for the elastic cross section. 

To elucidate their equivalence, we first note that the pole fac¬ 
torization 0 relies on the power series expansion of the Jost 
function around its zero (at k - p). Recall that we assume a 
short range potential, such that T(k) is analytic in the entire 
complex k plane, or at least inside a wide region containing the 
threshold k = 0. We thus write 

OO 

T{k) = Y J F nip-k) n , 

n -1 

where the coefficients F„ are complex numbers. The absence 
of the term for n — 0 in the sum above makes it explicit that 
r F{p) = 0, and allows for the factorization 'F(k) = (p - k)F(k). 


with all coefficients real. In the case of a NTR, Ao is vanishingly 
small, and to a good approximation we have 

„ _Ao 1 

K " B x a 

where a is the full scattering length, which agrees with Eq. ( ITTI i. 
As we discussed in Sec. EH the S matrix is simply written as 

_ Ao + iB\k + A2k 2 + --- 
’ “ l ••/; / • .\-A •••' 

When the resonance pole p - ik is very close to k - 0, this ex¬ 
pression truncated to second order becomes equivalent with the 
similar expression from the pole factorization approach, which 
has the advantage of using the exact value of the resonance pole 
(extracted as a fitting parameter from the computed scattering 
cross section, as it was done in Sec.0. However, for complex 
k, the very definition of the S matrix becomes problematic, as 
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it imposes strict requirements for the asymptotic behavior of 
the interaction potential. Hence, the direct approach based on 
Jost functions for real k has the advantage of being fully gen¬ 
eral. We also emphasize that the Jost function has a smooth k- 
dependence, and it can be easily interpolated on a rather coarse 
/.'-mesh, which is very convenient in practice; nevertheless, this 
advantage has remained largely overlooked, despite being rec¬ 
ognized more than four decades ago by Heller and Reinhardt 
& 

Regarding the coupled-channel case, we remark that Eq. (l57l > 
obtained in the Jost function approach represents a parametriza- 
tion which is identical to Eq. (IT71) from the pole factorization 
treatment, provided we simplify the latter to read k(r vcdcl (k) = 
Anq{k 2 + |/?| 2 )~ 2 . We can thus identify Dq/Di = \p \ 2 = q 2 + k 2 
in the denominator, and we regard q in the numerator as an 
overall measure of background reactivity (due to the couplings 
Vni = Vin between the entrance channel and all other channels). 

7. Qualitative interpretation based on the wave function am¬ 
plitude 

A very simple picture emerges from the treatment using the 
Jost function; indeed, as noted for the single-channel discus¬ 
sion in Sec. ED the amplitude of the physical wave function 
is modulated by the inverse of the Jost function J r ~ i {k), see 
Eq. (l28l >. This scaling property is especially useful in the low 
energy domain, i.e., for k ^ \a\~ l , when the regular solution 
in the short range (R F |a|) region is nearly ^-independent: 
<f>k(R) ~ <!>t)(R)- Hence, within the short-range region, Eq. (l28l) 
reads t/q-(R) ~ kkF~ l (Jc)<f>o(R), and the k dependence of the phys¬ 
ical solution is driven entirely by the scaling factor k r F ] (k). 

In the absence of NTRs, the Jost function can also be re¬ 
placed by its value at k — 0, 0) = Aq, and we obtain the 

well known linear scaling of the short-range amplitude of the 
physical wave function, i//k(R) ~ kAk x <po{R), characteristic to 
the Winger regime. However, when f F(0) = Aq is vanishingly 
small, it produces a resonance enhancement for if/k{R), and also 
the anomalous NTR behavior, 

WR) \ 2 ^ A 2 + B]k 2 lMR f - (58) 

We remark that this simple expression for the single-channel 
case already contains the main result, Eq. (l57l) . obtained rig¬ 
orously in Sec. 15.21 for the coupled-channel case. This can be 
justified qualitatively, if we rely on the following physical ar¬ 
gument; under the assumption of small couplings, the entrance- 
channel component (t/q) of the physical wave function follows 
the same scaling as in the single-channel case, see Eq. G3- 
This /.-dependence in the equation above will be imprinted (via 
the couplings V nm ) on all the other components i//„. Although 
the entrance-channel component obeys Eq. (l58l > only at short 
range, all other components will obey it for all R (including the 
long-range region). Thus, according to Eq. ( 1381 ) for m = 1, we 
extract the matrix elements S „ i from the asymptotic behavior 


of iJ/n(R), and we obtain for n 4- 1 


T~\Sm\ x\ifr n \ oc ———, 
kn A- + B 2 k 2 


where we used k - k\ for the entrance channel, and k„(k) x 
is constant. Finally, we obtain for the reaction cross sec¬ 
tion. 


kcr KdCl (n <- 1) ■ 


A 2 ■ 


■ B~k 2 ’ 


which is the same result obtained in Eq. (l57l) . with D o = A 2 and 
D 2 » B 2 . 


8. Conclusions 


In conclusion, we found that a near threshold resonance will 
affect dramatically the behavior of the ,v-wave cross section at 
low energy. In fact, the low energy domain will be divided in 
two distinct regimes: the Wigner regime with the well known 
scaling (cr react a- k 1 ) and the NTR regime with the anomalous 
scaling <r rcact oc ZT 3 . We derived simple analytical expressions 
for both elastic and reaction (inelastic) cross sections, which de¬ 
pend on the position of the pole of the S matrix, and found very 
good agreement with results obtained numerically for the full 
reactive scattering problem. We explained the anomalous NTR- 
regime behavior using two different approaches—one based on 
the pole factorization of the S matrix, and a more general treat¬ 
ment based on the low energy expansion of the Jost function. 

We remark that the Zr 3 NTR-regime behavior (for s-wave) 
is a general feature, and it can be a rather common occurrence 
in scattering problems; however, unless the resonance pole is 
very close to k = 0, the presence of the NTR can be somewhat 
inconspicuous, or it can be masked by higher partial wave con¬ 
tributions. Although previous work hinted at such anomalous 
resonant behavior—which was explored either by using mass¬ 
scaling inm, or external fields 149]—the new NTR-re gim e 
behavior shown in Eq. ( fl5l ) was discussed only recently [24]. 
These studies explored the effects of mass-scaling or magnetic 
field scanning at fixed scattering energy, e.g., see [50]. In our 
work, we revealed the effect of NTRs in two benchmark atom- 
diatom reactive scattering systems by mass-scaling (Hi+CI and 
Ht+F), but it also appears in any system with a zero-energy 
resonance, such as in photoassociation at ultralow temperatures 
1 5lLl52jl , ultracold collisions in general [53], or spin-relaxation 
in ultracold atomic samples 15411 . The modified (/T 3 ) behavior 
of the reaction (inelastic) cross section will play a role in the 
interpretation of experiments such as [fjj, and also for theo¬ 
ries based on Wigner’s threshold law developed to account for 
resonances in ultracold molecular systems fll7ll . 

This work was partially supported by the US Department of 
Energy, Office of Basic Energy Sciences (RC), and the Army 
Research Office Chemistry Division (f. S.). 
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